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Heating in nanoscale systems driven out of equilibrium is of fundamental importance, has ramifications for 
technological applications, and is a challenge to characterize experimentally. Prior experiments using 
nanoscale Junctions have largely focused on heating of ionic degrees of freedom, while heating of the 
electrons has been mostly neglected. We report measurements in atomic-scale Au break junctions, in which 
the bias-driven component of the current noise is used as a probe of the electronic distribution. At low biases 
(<150 mV) the noise is consistent with expectations of shot noise at a fixed electronic temperature. At 
higher biases, a nonlinear dependence of the noise power is observed. We consider candidate mechanisms 
for this increase, including flicker noise (due to ionic motion), heating of the bulk electrodes, 
nonequUibrium electron-phonon effects, and local heating of the electronic distribution impinging on the 
ballistic junction. We find that flicker noise and bulk heating are quantitatively unlikely to explain the 
observations. We discuss the implications of these observations for other nanoscale systems, and 
experimental tests to distinguish vibrational and electron interaction mechanisms for the enhanced noise. 

Nanoscale junctions are ideal tools for investigating the fundamental processes of heating and dissipation 
in many-body systems driven out of equilibrium, a situation fraught with complications. Figure 1 shows 
some of the relevant lengthscales to consider in this problem. A ballistic constriction (an atomic-scale 
junction') bridges between macroscopic source and drain electrodes. With no bias or temperature gradient 
applied, the electronic distributions in the source and drain and throughout the junction region are described 
locally as Fermi-Dirac distributions, characterized by some electrochemical potential (Fermi energy) and a spread 
in energies given by k^Tg, the equilibrium temperature of the system. Likewise, the local phonon populations 
(including optical phonons) are characterized by the same temperature, Tq. Finally, any local ionic excitations 
such as the ubiquitous tunneling two-level systems (TLS) undergo dynamics determined by their coupling to the 
equilibrium phonon and electronic distributions. 

When a voltage bias is applied and the system is allowed to reach steady-state, the situation is considerably 
more complicated, as shown in Figure 1. Far from the junction (region A), the widths in energy of the electronic 
and vibrational distributions remain described by k^To, presuming that heat is transported away into the bulk 
effectively. However, the Fermi energy of the source is raised over that of the drain by an amount eV. In the 
conventional treatment of such junctions^-*, the local electronic distribution functions within regions C and D are 
no longer simple Fermi-Dirac distributions. In the drain, for example, right-moving carriers that originated in the 
source (with energies ~ above the Fermi level of the drain) and transmitted through the ballistic region are 
present. However, it is generally assumed that those transmitted carriers are drawn from a Fermi-Dirac distri- 
bution with a temperature a Tq. We note that dynamical TLS and mobile atoms in region C that couple to 
electrons are sensitive to the full range of available electronic energies, and hence can respond'*'^ as if coupled to an 
effective temperature comparable to eV, generally much larger than To or T^. On a length scale set by the electron- 
electron inelastic mean free path (lee)- the electrons can exchange energy among themselves, evolving their 
distribution toward a Fermi-Dirac distribution characterized by a position-dependent, elevated electronic tem- 
perature > Tq. On a still longer distance scale, the electron-phonon inelastic mean free path (Le-ph)> the 
electrons and phonons can both reach approximately thermal distributions given by position-dependent tem- 
peratures Te = Tph > To- (The heirarchy of these length scales depends on relative inelastic scattering rates'". The 
assumption L^^ < Le-ph is reasonable in Au at room temperature at a bias of 0.1 V', and L^e is likely to be further 
reduced in a junction's confined geometry".) 

In general, it is difficult to define "temperature" for a system driven out of equilibrium and indeed several 
definitions can be advanced'. In the above, the electronic and vibrational distribution functions in regions A and B 
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Figure 1 | A nanoscale junction driven out of equilibrium. On length scales long compared to the mean free path for inelastic electron-phonon scattering 
(beyond region A) the electronic and ionic degrees of freedom are assumed to approach approximately thermal distributions parametrized by identical 
(position-dependent) temperatures (elevated above the equilibrium temperature, Tq, by an amount depending on the thermal path to "the 
environment"). In region B, the electronic distributions quasi-equilibrate into a local Fermi-Dirac thermal form, parametrized by a position-dependent 
electronic temperature, T„ greater than that of the ionic degrees of freedom. In region C, close to the ballistic junction relative to the mean free path for 
inelastic electron-electron scattering, the electronic distributions are non-thermal, resembling a linear combination of Fermi-Dirac distributions each 
with some electronic temperature Tj and energetically offset from each other by the applied bias, eV. Transmission through the ballistic region is through 
channels of transmittances t,-. Quantum corrections to the ballistic conductance can arise through quantum interference of trajectories (black and red) 
involving scattering off disorder in the diffusive electrodes within a coherence length of the ballistic region. The nonequilibrium electronic distribution in 
regions C and D can lead to nonequilibrium populations of local phonon modes. 



are well described as having the thermal functional form, with the 
effective temperature defined through that description. In regions C 
and D, the electronic and ionic populations are not thermally dis- 
tributed. However, it is still possible to parametrize such a nonther- 
mal distribution with a single effective temperature. For example, 
one can consider allowing local interactions with a large auxiliary 
"bath" (a thermal probe), and determining Teff as the temperature of 
the bath such that there is no net steady-state energy flux between the 
bath and the system'. Taking an experimental observable (bond 
strength; current noise; etc.) and comparing that observable with 
an empirical equilibrium temperature dependence or with a theor- 
etical model containing a temperature parameter is a long-estab- 
lished approach to defining an "effective" temperature, though one 
should be wary of over-interpreting the results. 

In the conventional treatment of these systems^'^, the electronic 
distributions impinging on the ballistic constriction (region D) are 
well approximated as Fermi-Dirac distributions with relative Fermi 
energies shifted by eV. There can be corrections to this due to back- 
scattering from disorder in the diffusive source and drain'" and 
these are modeled as corrections to the overall transmittance. In the 
standard picture, the effective electronic temperature of those FD 
distributions is established by energy redistribution among the elec- 
trons via inelastic electron-electron scattering on the scale of Lee fn 
region B, and by balancing the influx of energy from the bias-driven 
nonequilibrium distributions and the outflow of energy by thermal 
conduction of the electrons and eventually the phonons. The essen- 
tial work by Henny et al." shows that at cryogenic temperatures it is 
easily possible to get significantly elevated electron temperatures 
(that is, Te — To > To) varying over tens of nanometers. 

There have been a growing number of experiments to characterize 
the energy distributions of the electronic and ionic degrees of free- 
dom in such nanoscale structures. In the case of ionic degrees of 
freedom, vibrational occupancies can be inferred directly by highly 
local Raman spectroscopy measurements" Estimates of effective 
ionic temperatures may be obtained through measurements of rup- 
ture strength of atomic-scale junctions'""^". Ionic temperatures may 
also be determined by measuring the response of an auxiliary system 
designed to exchange energy with the ions, such as a resistive 



thermometer^' or a nanoscale thermocouple^^, though these typically 
indicate effective temperatures tens of nanometers away from the 
junction. 

It is difficult to access the electronic distributions directly, particu- 
larly the nonthermal distributions in regions C and D. These non- 
thermal electronic distributions have been inferred quantitatively in 
some elegant experiments in metals^''^* and carbon nanotubes^"''^'" via 
local tunneling measurements. Noise measurements can provide 
comparatively direct access to the electronic distribution functions 
impinging on the constriction (the Fermi-Dirac distributions that go 
into the nonequilibrium distribution in region D, which dominates 
the transport). In equilibrium, the relationship between fluctuations 
and dissipation leads to the Johnson-Nyquist contribution^''^", cur- 
rent noise that is white over a broad frequency range with a noise 
power proportional to the absolute (electronic) temperature and the 
electrical conductance, Sj = 4:k^Ti.G A^/Hz. Nyquist-Johnson noise is 
well established as a primary thermometer^'. Additional "excess" or 
"shot" noise can appear when a nanostructure is biased out of equi- 
librium, due to the discrete nature of the electronic charge. Previous 
shot noise experiments" in the diffusive limit, where a metal nanos- 
tructure is large compared to the mean free paths for both elastic and 
electron-electron inelastic scattering (Lee)> have found that the 
Nyquist-Johnson noise and shot noise clearly indicate elevated, 
bias-driven, position-dependent electronic temperatures T^{V) > 
To in mesoscopic wires at cryogenic temperatures. Those cryogenic 
experiments, accompanied by modeling, demonstrate that this heat- 
ing is a result of the comparatively large thermal resistances of the 
reservoirs (regions A, B). We must be mindful of this reservoir heat- 
ing in any such experiment. 

Beyond this quasi-equUibration of the electron subsystem at some 
effective Te that takes place in region B, it is an open question whether 
electron-electron interactions elevate the effective electronic temper- 
ature in the immediate vicinity of the junction (regions C and D). In 
other words, is the energetic width of the (approximate) FD distribu- 
tions impinging on the ballistic region D the same as that at the 
nearby edge of region B, or are there local processes that further 
broaden those distributions (corresponding to an effective local 
increase in that is not the reservoir heating described above)? 
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Theoretical predictions exist regarding local increases of the elec- 
tronic temperature' '" as a function of bias, but experimental probes 
have thus far been indirect'^'^". For instance, the bias evolution of a 
continuum background in anti- Stokes Raman scattering'^ ascribed 
to the electrons has implied effective electronic temperatures ele- 
vated by as much as several hundred Kelvin at biases of 0.4 V. 
Rupture strengths of molecule-containing junctions, accessing the 
electrons through their effects on the stability of ionic bonding, were 
reported^" to be consistent with modeling of electronic heating in a 
hydrodynamic treatment of the electronic fluid'^ '". 

For a junction with discrete quantum channels ; having transmit- 
tances t„ if the source and drain electronic distributions are thermal 
with a temperature T, the expected form of the excess noise is"": 



Si — Go 



4hTj2^^+2eV cothf-^) ^t,(1-t, 



2hT 



(1) 



where Go = 2e^ /his the quantum of conductance, including the spin 
degeneracy. At zero bias this simplifies to the Johnson-Nyquist noise, 
while at zero temperature this reduces to: 



S/ = 2eyGo^T,(l- 



(2) 



The finite-T expression in Eq. (1) has been confirmed at cryo- 
genic'"''^ and room temperature"' in mechanical break junctions. 
Often the excess noise is characterized by the Fano factor, 

F= 1,(1—1;)^/^^. T,, the ratio of the (zero temperature) 

bias-driven excess noise to the Schottky case" of electrons arriving 
in Poissonian fashion, S/ = 2 eVG. 

Inelastic processes coupling electrons and local vibrational modes 
can modify p"-**. In the limit where electron-phonon scattering can 
significantly cool the electrons, F ^ 0 as in macroscopic conduc- 
tors*'. However, depending on the particular transmittance of chan- 
nels, F may be either enhanced or suppressed once eV exceeds the 
level splitting of the local phonon mode, ha). These changes have 
been observed experimentally", though in multichannel junctions 
unraveling the details can be complicated'". Depending on whether 
the local phonon subsystem can equilibrate with bulk phonons 
quickly, or whether the local phonon population is pumped by the 
electrons, various dependences of F on the applied bias V are 
possible''^"'". 

In this work we report room temperature measurements of the 
bias dependence of the current noise in scanning tunneling micro- 
scope-style Au break junctions, through a lock-in method based on 
broad band RF detection. We move beyond previous work"' *^ with a 
much larger volume of data obtained at biases as much as three times 
higher, specifically to consider the regime examined in other experi- 
ments"'^" that report circumstantial evidence of local electronic heat- 
ing''" in atomic- and molecular-scale junctions. At comparatively 
low biases ( I V| < 150 mV) we find bias scaling consistent with excess 
noise as described in Eq. (1), with an effective electronic temperature 
equal to the environmental temperature. At higher biases, the mea- 
sured noise increases with bias more rapidly than expected from Eq. 
(1). A simple, conservative estimate shows that bulk heating of the 
electrodes, definitely relevant at cryogenic temperatures", is com- 
paratively negligible in these 300 K experiments. We consider flicker 
noise as a possible contributor to the measured effect and find that 
this is unlikely to be quantitatively compatible with the observations. 
We examine whether the increase in noise could be consistent with 
nonequilibrium electron-phonon effects and/or estimates of elec- 
tronic heating local to the ballistic region. 

Results 

Our experimental approach has been published previously"' *^ and is 
described in detail in Methods. A gold tip is moved repeatedly at a 



steady rate towards and away from an evaporated gold film, such that 
atomic-scale Au junctions are made and broken cyclically. We per- 
form electrical measurements during this process, cycling an applied 
voltage Vo across the junction (and a series resistor) and measuring 
the corresponding current. A lock-in approach simultaneously gives 
the change in integrated radio frequency noise, proportional to Si{V) 
— S/(0), where V is the voltage drop across the junction. The radio 
frequency band employed in the experimental analysis here is from 
250 to 600 MHz; data taken with a bandwidth from 400 to 800 MHz 
is qualitatively identical. To obtain the ensemble-averaged noise as 
a function of bias at a given conductance, conductance and noise 
histograms (500-1000 breaking traces) are acquired at many applied 
voltage levels. Figure 2 shows a representative conductance 
histogram, along with ensemble-averaged noise as a function of 
conductance acquired at several different bias voltages. The 
ensemble-averaged data are then binned by conductance and plotted 
as a function of V. In the picture of Eq. (1), for bias-independent 
electron temperature and a given junction configuration, plotting 
Si{V) — Si{V = 0) as a function of the scaled bias 
eV . eV 

should give a straight line with a 



4ksTG 



[2ksT 
slope that is 

Ti(l-T,-)/ 



coth 



Fano factor 
to find the 



2hT 

precisely the zero-temperature 
y Zi{\ — Xi) 1 2_^. T,- In our case we expect 

ensemble-averaged Fano factor, (^^^- — t;) /X^- ^7 

EN 
Tj, SO that the slope 

should be^^%,(l-T,)y(G/Go). 

Figure 3(a) shows low bias linear fits at several different conduc- 
tances away from the peaks in the conductance histograms. As 
reported previously"", small nonlinearities in these scaled plots were 
observed for conductances close to quantum suppression points, 
ascribed to extra bias-induced excess noise caused by flicker noise 



12000 



8000 



4000 



(a) 













Figure 2 | (a) An example of a typical conductance histogram as well as the 
noise measurement. The peaks in histogram indicate conductance 
quantization, coincident with the relative stability of particular junction 
geometries. In single valence electron materials like gold, this conductance 
quantization coincides with suppression of the shot noise, indicating that 
the chemically stable geometries include significant contributions from 
fully transmitting channels, (b) Excess noise as a function of conductance 
for a series of bias voltages (from bottom to top, 40, 100, 160, 200, 260, and 
320 mV). The vertical lines indicate the particular conductances (purple = 
0.8 Go, blue = 3 Go, red = 5 Go) for which noise as a function of bias is 
plotted in Fig. 3. 
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(and possibly some electronic heating of the type discussed below) 
while shot noise is minimized at suppressions. 

After further increasing bias to higher levels, as shown in Fig. 3(b), 
significant nonlinearities as a function of scaled bias are present at all 
the conductance values, well above the extrapolated linear depend- 
ence found at low bias. Here different markers represent different 
independent data sets while each color indicates a particular conduc- 
tance value. The solid lines are the low bias linear fits. In all the data 
sets, the nonlinearity is concave upward and independent of the 
order in which the voltages were applied (high to low; low to high; 
or interleaved). This indicates that the nonlinearity is an intrinsic 
function of the bias, not a result of irreversible changes to the junc- 
tions during the sequence of histogram acquisition at different 
applied biases. The rather similar Fano factors in the several Go range 
are consistent with our earlier measurements"", and may imply the 
involvement of more quantum channels, suggesting an approach to 
noise properties similar to those in the diffusive regime'"'"^'. 

Discussion 

At a given conductance, at comparatively low bias, the measured 
noise scales with bias as expected from Eq. (1). At comparatively 
high bias, the measured noise increases more rapidly than this 




0 0.5 1 1.5 2 2.5 3 



4KTG[eV/2KTcoth(eV/2KT)-1] (10 "a^/Hz) 

Figure 3 | Excess noise power vs the scaled bias at (a) comparatively low 
bias, and (b) comparatively liigh bias. The shaded region in (b) is the 
domain shown in (a). Colors represent data at different conductances 
(purple = 0.8 Go, blue = 3 Gq, red = 5 Go), as in Fig. 2. The solid lines are 
corresponding linear fits. Ensemble-average Fano factors at 0.8 Go, 3 Go, 
and 5 Gq are 0.58, 0.35, and 0.32, respectively, (a) At low biases away from 
the noise suppressions, the excess noise power is quite linear as a function 
of the scaled bias, consistent with shot noise as in Eq. (1). (b) At higher 
biases, significant nonlinearities are present in excess of extrapolations of 
the low bias linear dependence on the scaled coordinate. Note that 
independent data sets (acquired after tip cleaning and annealing) plotted 
in different markers give consistent results for these nonlinearities. The 
fiUed circle symbols correspond to the data sets in (a) and in Fig. 2(b). 



expectation. We must consider possible explanations for this trend 
in the noise at high bias. Possible contributors include: flicker noise; 
significant heating of the electronic reservoirs (region B in Fig. 1); 
electron-phonon inelastic corrections to the noise; and local elec- 
tronic heating (in region C of Fig. 1). 

Flicker noise is produced by fluctuations in the device resistance 
due to scattering of the electrons by dynamical defects^^"^''. At low 
bias, such that the dynamical defects remain effectively at Tq, flicker 
noise manifests as a voltage noise power that scales quadraticaUy in 
the applied dc bias. At higher biases, because of the availablity of 
carriers as much as above the local equilibrium Fermi level, the 
dynamical defects can have effective temperatures (of their ionic 
degrees of freedom) that are considerably elevated'*. The result of this 
ionic heating is increased flicker noise as more defects are able to 
participate, and a bias dependence of the voltage noise that is super- 
quadratic^^''"''. The operative question is whether the increased noise 
we see at high bias is indicative of this kind of enhanced flicker noise 
contribution. 

We assess the role of flicker noise in multiple ways: considering the 
overall contribution of flicker noise to the total noise signal; exam- 
ining the scaling of noise response with RF frequency range; and 
examining the expected dependence of noise on bias and junction 
conductance. The noise measurements in this work are broadband 
and in the rf regime, a frequency range (250 MHz to 600 MHz) 
considerably higher than that in which flicker noise is typically mea- 
sured. The presence of strong noise suppression at 1 Gq at the highest 
biases (see, e.g.. Fig. 2) establishes that the bulk of the measured noise 
likely results from shot noise of the type described at low bias by Eq. 
(1). We have also performed noise measurements with a different RF 
filter set (limiting the band to between 400 MHz and 800 MHz) at a 
variety of biases. While environmental backgrounds proved more 
annoying over that higher frequency band, the relative magnitude 
of the noise suppression near values of quantized conductance is 
essentially unchanged. This again is consistent with the significant 
majority of the total measured noise being shot noise, since one 
would expect overall reduced flicker noise at the higher frequencies. 

We look at the bias dependence of the noise. In previous measure- 
ments^*", when looking at 1 Gq, where the shot noise contribution is 
maximally (but not completely) suppressed, we found residual non- 
linearity with bias in the ensemble-averaged noise power at low bias 
that was roughly consistent with a flicker noise contribution. We 
repeat such an analysis here. Note that the superlinear behavior in 
Fig. 3(b) is as a function of scaled bias, not just V. Figure 4 shows the 

3 I • < 1 • • • 1 




Bias Voltage (V) 

Figure 4 | Excess noise at 3 Gq converted into voltage noise power Sv^as a 
function of bias. The solid blue line is the best quadratic fit to the 
experimental data, while the dotted blue line is a systematic overestimate 
(see text) of the contribution of flicker noise at 1 Go. Any flicker noise 
contribution at 3 Go is expected to be even smaller. Thus, flicker noise is 
unlikely to be a compatible explanation for the observed increase in noise 
as a function of V. 
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measured noise (converted to voltage noise, Sy) as a function 
of bias across the junction, V, for G = 3 Go- The data are well 
described as a linear + quadratic term. There is no clear evidence 
of any superquadratic dependence of the type expected for high bias 
flicker noise. 

Finally, we consider the expected dependence of noise with the 
junction conductance. We can be overly conservative and ascribe all 
of the nonlinearity in Sy as a function of V to flicker noise. This 
knowingly overestimates the magnitude of flicker noise at 1 Gq, since 
the true finite-temperature expression for shot noise (Eq. (1)) 
has some nonlinearity due to the coth(eV/2/cBT) factor. Previous 
experimental work at lower frequencies (<400 kHz)^' used litho- 
graphically defined mechanical break junctions at cryogenic tem- 
peratures to examine flicker noise down to conductances 
comparable to Go in individual junction configurations. In those 
experiments, for a fixed bandwidth, the investigators found that 
the ensemble-averaged flicker noise scales like Sy/V^ oc G"'"'' as 
junctions approach the ballistic regime, though there is variation 
within the ensemble. In Fig. 4, we compare the bias-dependent noise 
at 3 Go with that at 1 Gq (dashed line). Contrary to the ensemble- 
average scaling expectation for flicker noise, we find that the 
observed nonlinearity at 3 Go is nearly two times larger than the 
overestimate of the 1 Gq flicker noise. While large variations are 
possible between the noise properties of individual junctions^', this 
trend (greater nonlinearity vs. V and noise magnitude at higher 
conductance than at lower conductance) is robust throughout our 
ensemble-averaged measurements, in both RF frequency bands used, 
and goes against expectations for flicker noise from previous stud- 
ies^^. Thus, based on limits to the overall contribution of flicker noise 
to the total measured signal, and the systematic noise response as a 
function of frequency band, bias, and junction conductance, we find 
it quantitatively unlikely that the increase in noise observed at high 
bias in Fig. 3(b) is a result of flicker noise. 

A second possible explanation for the observed data that must be 
considered is heating of the "reservoir" electrodes due to current 
flow. In prior experiments by Henny et al." considering shot noise 
at cryogenic temperatures in a diffusive wire connected to thick fUm 
reservoirs, the heating of those reservoirs was definitely non-neg- 
ligible. In the present experiment, the situation is considerably dif- 
ferent: an atomic-scale junction in what is usually considered the 
ballistic regime, at room temperature, with most of the dissipation 
happening in effectively bulk reservoirs. The Au tip has a measured 
electrical resistance of less than 1 Q, and the sheet resistance of the 
Au film is similarly around 2 Q. The question is whether power 
dissipation near the junction could, through the finite thermal res- 
istance of the electrodes, lead to significant heating of the electrodes 
(both electrons and phonons) near the junction. To be very conser- 
vative, assume that only the electronic component of the thermal 
conductivity is relevant. Furthermore, assume that 100% of the dis- 
sipated power is deposited in either one of the electrodes proximal to 
the junction. For 300 mV bias across a 3 kfi junction, the power P = 
3 X 10"^ W, and this must correspond to (T^^^ — T^)LGres/ 2, where 
r^es is the local temperature of the reservoir electrode near the junc- 
tion. To is the environmental temperature, L = 2.44 X 10"" W fl/K^ 
is the Lorenz number, and Gres is the conductance of the reservoir, in 
this case one Siemens. For To = 300 K, we find T„s = 304 K. In the 
real situation, when roughly half of the power will be dissipated in the 
other electrode^^ and phonons wiU contribute significantly to 
the thermal path in both electrodes, the actual local temperature rise 
of the bulk electrodes (region B in Fig. 1) is expected to be even 
smaller. Note that this simple estimate based on the electronic ther- 
mal conduction is independent of the detailed electrode geometry. 
Therefore, we conclude that this kind of reservoir heating, while 
definitely important in the cryogenic experiments of Henny et aU^, 
is unlikely to be a significant contributor to the high bias noise 
increase in the present work. 



At the biases considered in this work, and at room temperature, it 
is important to consider the possible effects of electron-vibrational 
coupling on the measured noise. As mentioned previously, there are 
many theoretical predictions describing electronic transport in small, 
ballistic junctions coupled to a local vibrational mode'"**''''. In an 
individual junction one expects a modification in the Fano factor 
of the shot noise as the bias exceeds the characteristic energy of the 
vibrational mode. The change in P depends sensitively on the par- 
ticular channel transmittance'^ '"'''"'''*, as well as the coupling of the 
channels to the phonon degree of freedom. 

A detailed comparison with particular theoretical descriptions is 
difficult: Our data is in the ensemble average, generally with multiple 
contributing channels; if one expects the Au optical phonon to be a 
dominant contribution, our data is acquired at the limit kgT> hw so 
that finite thermal phonon population is possible. A detailed consid- 
eration of ensemble averaging over channel distributions would have 
to match the observed increase in enhanced noise at higher conduc- 
tances. However, the observed power law dependence of the excess 
noise, ~ V^, experimentally constrains possible phonon-based 
mechanisms for the increase in noise at high bias. Descriptions that 
assume strong pumping of the local phonon population by the elec- 
trons predict voltage dependences with V and V* contributions'" ''*. 
For vibrational effects to be responsible for the enhanced noise, the 
ensemble-averaged system must be in the "equilibrated phonon" 
limit, such that the local phonon population N remains bias inde- 
pendent'*''. This suggests an experimental test for this mechanism. By 
performing such measurements as substrate temperatures are 
decreased down to the cryogenic regime, one would expect decreases 
in the local phonon population, and therefore decreases in such non- 
linear corrections to the noise. 

Given that temperature gradients within the bulk of the electrodes 
are small, an additional possible explanation for the observed trends 
in the measured noise would be some amount of effective electronic 
heating local to the junction (that is, electronic distribution functions 
impinging locally on the ballistic region (region D) that evolve non- 
trivially with bias for reasons other than bulk heating of the reservoirs 
(region B)). For example, the single-particle scattering formalism 
neglects electronic heating due to the electronic viscosity, inherently 
an effect of electron-electron interactions'*''". Before addressing 
this specific proposal, we attempt to assess whether the observed 
trend could be consistent with some broadened electronic distri- 
bution function (the idea that describing the components of 
the nonequilibrium distribution function in region C of Fig. 1 is 
bias-dependent). To do this, we parametrize the true nonthermal 
carrier distribution at the junction location through an effective 
electronic temperature, r(V), that depends on the bias, appearing 
in Eq. (1). Obviously, a general nonequilibrium distribution does not 
have to be parametrized by a single effective temperature; however, 
some kind of parametrization is necessary to perform an analysis, 
and this is closely related to the approach of Henny et al. when 
considering "traditional" electronic heating. We employ the 
Landauer transmittance expression and the resulting Eq. (1), but 
assume that T = T{V) = Tq + dT{V), where To = r(y = 0) is the 
ambient temperature, ~ 300 K. To further simplify the analysis, we 
leverage the fact that eV':^2ksTo in the high bias regime of interest. 

eV 

We can then move forward approximating coth ^ —> 1, which we 

will check for self-consistency (to make sure that this continues to 
hold for the inferred T{V)). Thus Eq. (1) for a single junction may be 
reduced to 



S/(y) = Go 



4fcB ( To -h <5 r) J] -H 2ey £ T, ( 1 - T,) 



(3) 



And at zero bias, we already know that 
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Figure 5 | A toy model calculation of transport properties. We consider 
different models of how particular transport channels with transmittances 
T; contribute to the conductance, (a) Each channel contributes in turn as 
conductance is increased, with no channel mixing, (b) An alternative 
model, with considerable channel mixing. In both panels, the dotted lines 
indicate calculated Fano factors, the solid blue lines denote > x'f, and 
the solid red line shows 2_^. '^i > all as a function of G. 



s,(y=o) = Go 



4fcBro^Tf+4fcBro^T,(l- 



(4) 



which is purely Johnson-Nyquist noise. Solving for 8T{V) 
'Si{V)-Si{V = 0) 



ST{V) = 



4hG 



eV 
2ku 



-To ]F 



r. 



N 7 



where S/( V) — Sj{V = 0) is exactly the excess noise we measure. The 
ensemble-averaged Fano factor, F, can be determined from the low 
bias data described above. The only missing information is the ratio 

EN 
i 



N 2 



For our ensemble-averaged measurements, by specifying G 



we then need to find the ensemble average 



GoE 



N 2 



WhUe in a 



single junction with Fano factor Fg, one can write the parameter of 
interest as 1/(1 — Pq)) note that in general one cannot assume that 

G 



GoE: 



= 1/(1 — F), where F is the ensemble average. 



We can consider models of channel mixing to illuminate possible 
values, as shown in Fig. 5. The dotted line represents the Fano factor. 
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Figure 6 | Within the model of a locally elevated electronic temperature 
at the ballistic junction, the inferred electronic temperature from Eq. (5) 
vs bias voltage. For data sets at 13 different conductances 3, 3.25, 3.5, 3.75, 
4, 4.27, 4.5, 4.75, 5, 5.25, 5.51, 5.73, 6 Go, the colors/markers, respectively, 
are: blue point, red point, green point, purple point, khaki point, blue star, 
red star, green star, purple star, khaki star, blue open square, red open 
square, green open square. (The apparent decrease seen in the shaded 
region is not physical and results from the failure at low bias of the 
approximation co^{eVI2k^'V{Y)) ~ 1.) At high bias, inferred electronic 
temperatures increase monotonicaUy with V, approximately independent 
of G. Error bars are dominated by the statistical uncertainty in the inferred 
Fano factor of the low-bias data. Dashed line is a fit to the simple model 
expression that assumes a local quasithermal electronic distribution, 
thermal transport by the electrons, and a fraction a of the power W 
dissipated locally in the junction, &T = (^T^ + 2ctV^ /lY^^ — To, with a = 
0.02. 



N y^'' T 

the solid blue line . t^, and the solid red line — j^—^ . Figure 5(a) 

Ei '^i 

represents one limiting case, no channel mixture at all, with each 
channel turning on linearly sequence as the conductance is increased. 
Figure 5(b) assumes much stronger mixing of the channels, so strong 
that there would be little detectable shot noise suppression in this 
limit. In both limits, the ratio of interest approaches a constant (one) 
quickly after the conductance is increased above 1 Gq. As the diffus- 
ive limit is approached", F^ 1/3 because of partially open channels, 



. E, 

meaning 



E- 

■ 1 .5. In either case, it is safe to approximate — 



Er^? Ef 

as a constant on the order of 1 when G is much above 2 Go (though 
the precise value introduces a systematic uncertainty in what fol- 
lows). Within this approximation, we can evaluate Eq. (5), and the 
results are shown in Fig. 6. 

The effective electronic temperature elevations in Fig. 6 are calcu- 
lated from the data sets used to produce Fig. 3(b). In this figure, 13 
conductance values are chosen, representative of conductances 

eV 



higher than Gq. The approximation coth 



2knT 



> 1 is only valid out- 



side the shaded region. At a bias of several hundred mV, within this 
model an effective electronic temperature elevation as high as 350 K 
is inferred. The error bars are dominated by the uncertainty in F, as 
inferred from of the linear fits in Fig. 3(a). These local electronic 
temperature elevations are the same order as those inferred through 
completely independent experiments employing indirect optical 
methods" and bond rupture measurements^". 

We use a deliberately simplistic model to check the reasonableness 
of these (57 values. The inferred effective temperature elevations are 
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Figure 7 | Schematic of rf noise measurement scheme. 



approximately independent of conductance, which is consistent with 
a local electronic heating picture assuming a balance between local 
power generation and heat current carried by the electrons'". In this 
local heating picture, some fraction a of the total Joule heating is 
dissipated local to the junction, Pi„ = alV = aV^G. (In the standard 
approach, a would be zero; given the general success of the standard 
approach, the expectation is that a should be small compared to one.) 
Since the electron-phonon inelastic mean free path is much longer 
than the junction size, the electrons must carry away the heat gen- 
erated. In the limit of local quasithermal energy distributions (so that 
a local electronic temperature is well defined) for the outgoing car- 
riers, the electronic thermal conductance is LTG. Balancing the Joule 
heating with the thermal current Pout= {T^ — T^^LG/l leads to a 

predicted local temperature of (T^ -|-2ay^/L)^^^. The observed data 
are roughly consistent with this for a ~ 0.02, not bad considering the 
simplicity of the model, which assumes the Wiedemann-Franz law 
and no dependence of the t, on bias. The inferred 5T seems to 
increase more rapidly than this dependence at higher bias; this sim- 
ple model may be inadequate, assuming for example that the trans- 
mittances t, are independent of bias (generally not a poor 
assumption for Au over this energy range^^'^**). 

We see that the observed noise trend could be explained by a 
relatively modest amount of electronic heating local to the junction 
(dissipating 2% of the total power in region C of Fig. 1, though in 
general one would expect a. to depend on the conductance). We now 
consider whether there is a physical mechanism that could contrib- 
ute this amount of local electronic heating. Indeed, similar to the 
local, nonequUibrium ionic heating due to electron-phonon interac- 
tions, the internal Coulomb friction of the electron liquid is predicted 
to create local electron heating'". Due to the large current densities in 
nanojunctions compared to their bulk counterparts the electron- 
electron scattering rate increases locally in the junction'. The 
underlying Fermi sea thus "heats up" locally due to this increased 
scattering, via production of electron-hole pairs whose energy needs 
to be dissipated away from the junction. This would manifest itself in 
the noise just as an elevated electronic temperature does so in the 
regime of Henny et al". The local effect was first computed in Au 
nanojunctions'" assuming no dependence of the electronic viscosity 
on the junction geometry, namely on the electron transmission. 
Using the viscosity of that work, namely 10"' Pa-sec for the electron 
liquid of density comparable to gold (r^ = 3), we obtain an estimate of 
(5 T that is about two orders of magnitude smaller than what is shown 
in Fig. 6. The drawback of neglecting the dependence of viscosity on 
junction geometry was however understood and corrected in Roy 
et aV, where the zero frequency viscosity f;(0) was evaluated as the 
product of the high-frequency shear modulus and a junction- 
specific momentum relaxation time t. The actual scattering potential 
used in Ref 8 to compute t was chosen precisely to model Au 



quantum point contacts. By referring to Fig. 4 of Ref 8 we see that 
r\{Q) ~ 10"'' Pa sec for the transmissions typical of these point con- 
tacts. Since the zero-frequency viscosity enters as the square root of 
its value in the effective electronic temperature (see Ref 30) we obtain 
a value of a = 0.05, a factor 2.5 higher than the one fitted in Fig. 6, a 
very reasonable agreement considering that the actual geometry of 
the junction is not known. This quantitative consistency between the 
modeled experiment and this theoretical treatment suggests that the 
electronic viscosity should definitely be considered as a potential 
contributor to the observed enhanced noise. Within this theoretical 
picture the effective electronic heating depends strongly on the trans- 
mittance, implying that such heating would be much less important 
in lower conductance structures such as molecular junctions, where 
it might be detected through inelastic electron tunneling spectro- 
scopy linewidths". We further note that this electron-electron 
interaction effect should be comparatively independent of the back- 
ground equilibrium temperature. Therefore, measurements of the 
scaling of noise at high bias at lower substrate temperatures should 
allow the discrimination between electron- vibrational enhancement 
of the noise discussed above and this local electronic heating 
mechanism. 

In conclusion, we have conducted excess noise measurements in 
scanning tunneling microscope-style Au break junctions at room 
temperature and observed a nonlinear bias dependence of the noise 
exceeding the expectations for shot noise in the Landauer formalism 
at fixed electron temperature. The magnitude of this extra noise is 
inconsistent with flicker noise or heating of the bulk electrodes as 
possible confounding effects. NonequUibrium electron-vibrational 
effects could be relevant, but the observed power law of the noise 
increase suggests that any participating local phonon modes are not 
pumped significantly by the electrons. Interpreting the increased 
noise in terms of a model of local electronic heating, we find that 
the additional noise implies an effective electron temperature 
(defined through the width of the FD distributions impinging on 
the ballistic junction) elevated by as much as 300 K at 0.3 V bias. 
A theoretical treatment accounting for the effects of electronic vis- 
cosity is shown to be roughly consistent with such a value. It is 
important to note again that this noise method is only one way to 
parametrize the local nonequUibrium electronic distribution by an 
effective temperature. The temperature determined this way may not 
be identical to a local, nonequUibrium temperature defined'^ through 
the lack of net average transport of energy between a reservoir and 
the nonequUibrium system. Additional experiments looking at elec- 
tronic heating as a function of environmental temperature and trans- 
mittance should further constrain theoretical models of this 
nonequUibrium phenomenon, and make it possible to distinguish 
between electron-vibrational and electron- electron contributions to 
the enhanced noise. 
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Methods 

Measurement technique. A schematic of the measurement system is shown in Fig. 7. 
A function generator at 10 kHz cycles the applied bias between zero and a particular 
level Vo- This low frequency bias passes through a bias-tee and is applied to the tip; the 
film is connected, via the DC port of another bias-tee, through a resistance standard 
(as a current limiter) set to i^o — 2 kf2, and then to a current preamplifier. The current 
preamplifier output is measured directly by a high speed data acquisition (DAQ) card 
as well as feeding into a lock-in amplifier synced to the 10 kHz bias. Knowing Rq and 
the current, we calculate the true voltage drop across the junction, V. The high- 
frequency port of one bias-tee is terminated at 50 Ohms, while the other bias-tee's 
high frequency port is connected to an amplifier chain, a filter set to define the 
bandwidth {either 250 MHz to 600 MHz, or 400 MHz to 800 MHz), and then to a 
logarithmic power meter. The power meter output is read directly by the DAQ and is 
also directed into a second {"noise") lock-in amplifier synced to the low frequency 
bias. This second lock-in measures the component of the radio frequency noise 
(integrated over the gain -bandwidth product of the system including the amplifier 
chain) due to the applied "DC" bias. The noise lock-in signal (measuring the mean 
square voltage signal on the power detector) includes a background (independent of 
junction conductance and bias conditions) due to amplifier noise that may be 
subtracted away from the ensemble -averaged power measurements. Measurements 
of the junction's RF reflectance over the bandwidth as a function of conductance are 
essentially featureless between 0. 1 Go and 8 Gq, showing that the efficiency of the out- 
coupling of RF noise from the junction varies little over that conductance range. 
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